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Abstract 

Domain wall fermions are denned on a lattice with an extra direction the size of 
which controls the chiral properties of the theory. When gauge fields are coupled 
to domain wall fermions the extra direction is treated as an internal flavor space. 
Here it is found that this is not the case for scalar fields. Instead, the interaction 
takes place only along the link that connects the boundaries of the extra direction. 
This reveals a richness in the way different spin particles are coupled to domain wall 
fermions. As an application, 4-Fermi models are studied using large N techniques 
and the results are supported by numerical simulations with N=2. It is found that 
the chiral properties of domain wall fermions in these models are good across a large 
range of couplings and that a phase with parity-flavor broken symmetry can develop 
for negative bare masses if the number of sites along the extra direction is finite. 
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1 Introduction 



Domain wall fermions JT| provide an alternative to standard lattice fermions and have 
already been used to formulate lattice gauge theories (for reviews see and references 
therein). Here domain wall fermions are used to study lattice theories of fermions and 
scalars. Such theories can be studied analytically using large N methods and applications 
of these theories requiring enhanced control over the chiral symmetries may be of interest. 
For example, such applications may involve 4-Fermi models or models with Higgs fields 
interacting via Yukawa couplings. 

Domain wall fermions (DWF) are defined on the sites of a space-time lattice with one 
extra direction. The dimension of this lattice will be denoted as d + 1. In the method 
of || the boundary conditions along the extra direction are free and as a result light 
fermion surface states develop on the boundary. The plus chirality components of the 
Dirac spinor are exponentially localized on one wall while the minus ones are on the other. 
If the extra direction has L s sites then the two chiralities overlap only by an amount that 
is exponentially small in L s . Therefore, the extra parameter L s controls the amount by 
which the regulator breaks the chiral symmetry at any lattice spacing. As a result, and 
contrary to staggered or Wilson fermions, the chiral and continuum limits are separated. 
Furthermore, in numerical applications the cost for recovery of the chiral symmetry is only 
linear in L s . These properties make DWF an attractive regulator in problems where good 
control of the chiral symmetries is needed. For example, DWF have been used in studies 
of the finite temperature QCD phase transition jj], (51 ^ [7[ 

One would have expected that the formulation of theories with DWF coupled to scalar 
fields should closely follow that of DWF coupled to gauge fields. However, this is not the 
case and as a result an interesting difference between the way the scalar and gauge fields are 
coupled to DWF emerges. The formulation and a discussion of this difference is presented 
in section 0. Two 4-fermi models, one with a Z2 x Z2 and one with an 577(2) x SU(2) chiral 
symmetry in three dimensions, are studied in section || using large N techniques. In section 
[| these models are studied using Hybrid Monte Carlo simulations with N = 2. Conclusions 
are presented in section ^. 

Using DWF in 4-Fermi models was also suggested in 11. For a recent implementation 



of overlap fermions to the gauged Gross- Neveu model see [10] . 



2 The fermion— scalar interaction 

The free DWF action in the formulation of || is: 

S = - ^(x,s)D(x,s;x r ,s')^(x',s r ) (1) 

x,x',s,s' 

with Dirac operator 

D(x, s; x', s') = 6(s - s')fl(x, x') + p-(s, s')6(x - x') + E ± (s, s'; x)6{x - x') (2) 

where 

1 d 

fl(x, x') = - J2 [(1 + + ft - x') + (1 - llx )6{x' + pL-x)\ + (m - d)6{x - x') (3) 

1 M =l 
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P R 6(1 - s') - 6(s - s') 
p L (s, s') = \ P R 5(s + 1 - s') + P L 6(s -1-s') 
P L S(L S -2-s')- 6{s - s') 



s = 

6(s-s') 0<s<L s -l (4) 
s = L. - 1 



E ± (s, s') = -m f [Q R (s, s') + Q L (s, s')} 



R,L 



1 ±75 



(5) 
(6) 



Q R (s, s') = P R 5(L S - 1 - s)6(0 - s'), Q L (s, s') = P L 5(0 - s)5(L s -1-s') (7) 



In the above equations mo is a five-dimensional mass representing the "height" of the 
domain wall. Among other things, the parameter m controls the number of light species 
For example, in four dimensions the theory is symmetric around mo = 5 and for 
< m < 2 corresponds to a theory with one light species, for 2 < m < 4 to a theory 
with four and for 4 < m < 6 to a theory with six | IT| . The parameter rrif explicitly mixes 



the two chiralities and as a result it controls the bare fermion mass of the <i-dimensional 
effective theory in a linear fashion. The DWF Dirac operator satisfies the identity 0: 



lh RD Flh R = D\ 



(8) 



with R the reflection operator along the fifth direction. As a result, the two species deter- 
minant is real and non-negative 



det Dp 1 



det D\D F 



(9) 



From a calculation of the propagator or from the lowest eigenvalue of the fermion matrix 



the fermion mass for the one flavor free theory is [12 



m e fj = m (2 — m ) mj + (1 — m 



< m < 2 



(10) 



From the free action one can see that there is more than one way one could couple a 
gauge field to DWF in a gauge invariant way. The most straightforward choice of intro- 
ducing a gauge field is to also define it in d + 1 dimensions and couple it in the standard 
gauge invariant way. However, it turns out that this choice does not lead to the correct 
theory in d dimensions. The correct way to couple the gauge field was arrived at by think- 
ing of the extra dimension as yet another internal flavor space [|T], |TB|. Then one has to 
introduce the gauge field only in the d dimensional space and couple it the same way to all 
fermion "flavors" . Obviously, this also results to a manifestly gauge invariant formulation. 
In particular, the Dirac operator is as in eq. but with p defined as: 



MA,,') = l£ Ui+^U^Six + f^-x^ + il-^Ulix^x' + fl-x) 

A n=i 

+ (m — d)6(x — x) 



ir 
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Therefore, although one can think of the extra direction as a new space-time dimension 
it is more natural to think of it as an internal flavor space. From this point of view, if one 
wishes to formulate a theory that in <i-dimensions would be a theory of a scalar field coupled 
to a fermion bilinear, one would couple the scalar field to all internal fermion flavors the 
same way, as was done for the case of a gauge field. However, this choice will not result 
to the desired theory in d dimensions. To see this consider the simple example of a target 
theory of fermions coupled to a real valued scalar field a in d dimensions with continuum 
action: 

S = J ll^d^q + rriqq + aqq + (3a 2 (12) 

If m = this model has a Z 2 x Z 2 chiral symmetry. It is clear that a large N or mean 
field analysis of this model will indicate that the propagator mass would be m+ < a > 
where < a > is the vacuum expectation value of the a field. Clearly a lattice formulation 
of this model using DWF should at least reproduce this simple result. However, if the a 
field is coupled to all fermion degrees of freedom the same way, as is the mass m , i.e. 

similar analysis will result to a propagator mass as in eq. [It], but 
with mo replaced by mo+ < o >. This is obviously the wrong result. Even if one was 
to accept this exponential behavior by redefining the a field this theory would have the 
peculiar property of changing the number of light species depending on the value of a. 
Clearly, in order to get the correct result one has to couple the a field in the same fashion 
as the rrif mass and not as the d+ 1 mass mo. Then the Dirac operator is as in eq. but 
with E ± : 

E ± (s, s'; x) = -[m f + a{x)\ [Q R (s, s') + Q L (s, s')] (13) 

In the resulting action the a field couples with the fermion fields only across the bound- 
ary links of the d + 1 dimension. That this is the correct action can be seen by observing 
that the interaction term in the action can be rewritten as: 

L s - l)a(x)P R V(x, 0) + V(x, 0)a(x)P L V(x, L s - 1) = q(x)a(x)q(x) (14) 

with: 

q(x) = P R *{x,0)+P L V(x,L s -l) 

q(x) = V(x,L s -l)P R + *(x,0)P L (15) 

The fields q, q correspond to the light d- dimensional fields [0. Since these are the only 
propagating fields in d dimensions and since, in terms of these fields the correct global 
symmetry is in place, the above action should describe the desired target theory. 

Another way to see that this is the correct action is to consider adding such an inter- 



action to Neuberger fermions fl4| . In that formulation the fermion fields are standard d 
dimensional fields and the interaction can be added unambiguously. Then if the discrete 
finite flavor version of that formulation is used one can directly recover the DWF action. 
Since the a field is added as a mass term, one can see by simple inspection that the above 
action is recovered. Furthermore, the connection can also be made exact by using the 



results of [[HJ where the (^-dimensional action arrived at by integrating all heavy fermion 
degrees of freedom is expressed solely in terms of the q, q fields. 
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It should also be noted that scalar fields were used in a different context [|l], [16| in an 
attempt to formulate chiral gauge theories using DWF. There, in order to maintain gauge 
invariance, the scalar fields were defined on a "slice" in the extra direction. 



With the interaction defined as in eq. 13 the transfer matrix of the model along the 



d + 1 direction [[13], |3| is independent of the a field and is therefore the free field transfer 
matrix. The a field is defined only on the boundary and therefore only affects the boundary 
operator O of 0]. In the case of a gauge theory the extra fermion degrees of freedom are 
heavy but numerous as the L s — > oo limit is taken. For this reason it was shown that 
they need to be subtracted and the way to do that for each gauge field background is to 
divide out a bulk factor (T^j of the form f(U) Ls where f(U) is some function of the gauge 
field U determined by the transfer matrix (in the formulation of @ f(U) = A max detS). 
Here the same bulk factor is present but unlike the gauge theory it does not depend on 
the background field. Therefore it can be factored out from the a path integral and as a 
result it becomes an irrelevant factor that will cancel out with the same factor coming by 
dividing by the partition function whenever an expectation value is calculated. Therefore, 
unlike the gauge theory, it is unnecessary to explicitly subtract it out. 

The extension to more general interaction terms is straightforward. The interaction is 
written in terms of the q, q fields and using eq. O is transcribed for the \I/ fields. For 



example, a target theory with SU(2) x SU (2) chiral symmetry has Dirac operator as in eq. 
| but with E ± : 

E ± (s, s'; x) = -[mf + a(x) + ij 5 ir(x) • r] [Qr(s, s') + Ql(s, s')] (16) 

where w(x) is the pion field and r are the Pauli matrices. 

The above results, although contrary to naive expectations, are straightforward. How- 
ever, the following observations can be made: 

1) For the case of scalar fields it is not natural to interpret the extra DWF direction as 
an internal flavor space since the interaction takes place only at the boundary of the extra 
direction. 

2) One could have thought that the introduction of an extra dimension would have lent 
some freedom in the definition of the DWF interaction with the scalar and gauge fields. 
However, this appears not to be the case. The DWF formalism and symmetry requirements 
seem to have naturally forced the interaction terms to be of this specific form. 

3) A picture with some richness seems to develop with different spin fields coupled to 
domain wall fermions in different ways. 



3 Large N analysis 



The analysis in this section closely follows the analysis in |17j of the same model using 
Wilson fermions. Also see [IS, [19|, |2(], [21]] for studies of 4-Fermi models with staggered 
fermions. 

Consider a 4-Fermi model in the continuum with SU(2) x SU(2) chiral symmetry and 
Lagrangian: 



C = q(ifi- m)q + — 



{qqf + (qi^rqf 



(17) 
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In the above expression all indices have been suppressed. The fermionic field q is a flavor 
SU(2) doublet and a color SU(N) iV-column vector. The Lagrangian is diagonal in "color", 
in contrast with the full QCD Lagrangian which is diagonal in flavor, r = {ti,t 2 ,t 3 } are 
the three isospin Pauli matrices, f)= 7^(9^, and m is the bare quark mass (if m 7^ the chiral 
symmetry is explicitly broken). As is well known, a Lagrangian density that is quadratic 
in the fermionic fields can be obtained by introducing a scalar auxiliary field a and three 
pseudoscalar auxiliary fields 7r = {71*1, 7r 2 , 7r 3 } coupled to the fermion bilinear: 

C = qDq-nfp^a 2 + tt 2 ) 

D = if)— m — a — 275T • 7r . (18) 

Here rif is the number of flavors and f3\ = 2n 1 f . Gi ■ Transcribing this to a Euclidian lattice, 
using DWF with fermion- scalar interaction as described in section Q and an even number 
of "colors" N one obtains: 



Z = J [d^d^dadn}e~ 



N/2 

S = J2\WD¥+W +N/2 D^ i+N / 2 \+nMa 2 + 7r 2 ) (19) 



i=i 



where D is defined in eq. ^| with E 1 - as in eq. [16| and use of the property in eq. || has been 
made. 

Following the standard large iV analysis, the saddle point equations (SPE) are calculated 
by assuming uniform saddle fields and small fluctuations of the form: 

, . 5cr(x) , , 8tz(x) . . 

<t(x) = a s + -±J-, tt(x) = tt s + —±J- (20) 

The saddle point equations are: 

^ReTr^Qs + Q L )\ + 2n f ^a s = (21) 

^ReTr[Dj\Q R + Q L )ij 5 r] + 2n f 1 n s = (22) 



or in terms of the q, q fields of eq. [K: 

^Re < qq > s +2n f (3 1 a s = (23) 
^Re < qi^rq > +2n f (3 1 n s = (24) 

where ^ 

01 = i = 2^NGl (25) 

The L s = 00 and the finite L s cases need to be studied separately since each case has 
a different symmetry group. 

I. The L s = 00 case. 
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If mf = then the DWF Dirac operator has exact chiral symmetry. If there is sponta- 
neous symmetry breaking one is therefore free to choose the direction of breaking at will. 
The standard choice is a s ^ and iz s = 0. For this choice the second SPE eq. |22] is 
trivially satisfied since the non-trivial flavor part of D~ l is proportional to T7r s . This case 
is identical to the case of a simpler model without iv fields and rif — 1. Such a theory has 
Z2 x Z2 chiral symmetry. Therefore, the following analysis can be trivially extended for 
the Z 2 x Z 2 model. 



Using the free propagator results of |L2| after some algebra the first SPE eq. |21] for 
L s = 00 results in: 

_ £ y a = o (26) 

V p p 1 + (a s + m f ) 2 z(p) 2 

where 

f = MP,? (27) 
«(p) = 1 - be~ a (28) 

d 

6= i_ mo + ^(i_ C os(p M )) (29) 
n=x 

1 + b 2 +p 2 

cosh a = — , < a (30) 

(J 

The factor z(p) plays the role of "selecting" Brillouin zones. For a given range of mo the 
factor z(p) is nonzero at the orgins of only certain Brillouin zones. To see this observe that 
at the orgin of any Brillouin zone p ~ and therefore eq. ^ has solution exp(— a) = b if 
\b\ < 1 and exp(— a) — 1/6 if |6| > 1. Therefore, from eq. EHl one can see that z(p) = 1 — b 2 



if |6| < 1 and z(p) = if |6| > 1. For example, if < mo < 2 then z(p) is non zero only in 
the zone with momentum components around zero. The condition for nonzero z(p) at the 
orgins of a Brillouin zone: 

d 

\b\ < l^m -2< ^(l-cos(p M )) <m Q (31) 

At=l 

is the same as the condition for the existence of normalizable states |I[ . For zero momentum 
z{p = 0) = m (2 - m ). 

One can see that eq. ^ for m/ = can have two solutions: one with er s = corre- 
sponding to a chirally symmetric phase, and one with <j s 7^ corresponding to a phase 
with spontaneously broken chiral symmetry. The critical value of f3\ is obtained from eq. 



26] in the limit cr s — * and is given by: 



^ = l^-f (32) 

It is interesting to notice the rather strong dependence of f3i c on m . Similarly strong 
dependence of the critical coupling on mo was found for the QCD finite temperature phase 
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transition H. As an example, (3\ c is plotted versus mo for a 6 3 lattice with antiperiodic 
boundary conditions along the time direction in figure 1. In the next section such lattices 
will be used to compare the large N expressions with numerical simulations. A similar graph 
can be obtained for a four- dimensional lattice but it extends from m = to m = 10 and 
is symmetric around uiq = 5. 
II. The finite L s case. 

When L s is finite the DWF Dirac operator breaks chiral symmetry explicitly even for 
rrif = 0. Therefore a s ^ and the remaining symmetry is the 577(2) flavor symmetry. This 
symmetry can break spontaneously resulting to a nonzero ir s . If this happens then from 



eq. 24 one can see that < qi^rq >^ and the parity as well as flavor symmetries are 



broken. This is the Aoki phase [22]. This phase has also been observed for this model with 



Wilson fermions |23|, ^M. The existence of this phase for DWF was also suggested in p3 
and it may be present in QCD with DWF . 



First consider the phase with 7r s = 0. Again the second SPE eq. ^ is trivialy satisfied 
and the following analysis is also valid for the Z 2 x Z 2 model at finite L s . The first SPE 
eq. B2I is: 



2 

p 

+ e-^- 1 ) [(m f + a s )A m + b(B + A 1 + A 2 )} 

+ e - 2 < L ^[(m f + a s )A l + bA m ] } = (33) 

where Aq, Ax, A2, A m are functions of the momentum, rrif + a s and e~ aLs and are the same 
as in [|T^] but with rrif replaced by rrif + a s . This equation can be used to calculate a s as 
a function of the other parameters of the theory. 

It is interesting to see how chiral symmetry is restored as L s increases. As can be seen 
from eq. ^ the decay coefficient a is independent of (3i and it only depends on m and 
the momenta. In figure 2 the minimum and the maximum value of Re(a) (if b < a has 
imaginary part ±«7r) obtained for different momenta is plotted as a function of tuq in the 
range < tuq < 2. The "spikes" of the maximum value are the singularities that occur 
when 6 = 0. As can be deduced from eq. the peak of the minimum decay occurs at 
mo = 2 — y/2 = 0.586 and is — ln(2 — y2) = 0.535. The decay rates between the minimum 
and the maximum values have no gaps. As L s increases eventually the only chiral symmetry 
violations that remain will be controlled by the minimum decay rate unless the observable 
is dominated solely by terms with decay rate close to the maximum. Since the minimum 
decay rate is approximately constant across the full range of mo there is no valuable option 
of tuning m in order to achieve better characteristics except perhaps around m = 0.586. 
As an example, a s versus L s is plotted in figure 3 for mj = and various m at f3\ = 0.3 
which is above the transition for all the m values. The slopes obtained from the larger L s 
points are slightly faster than the minimum decay rates of figure 2 and the largest slope is 
for mo = 0.6. 

Next consider the phase with ir s 7^ 0. The first SPE is as in the 7r = case eq. 
3~3| . However, since the non-trivial flavor part of D~ x is proportional to T7r s one can 
now eliminate tz s from the second SPE and obtain a second non trivial equation. These 
two equations can be used to determine a s and iz s . The full form of these equations is 
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complicated and not particularly illuminating. However, if only the leading order in mj- 
and e~ aLs is kept these equations can be written in the form: 



0i 



V 



E 



z 

15 



+ 0(1) 



m t 



E P [e~ aLs 


( f 


V)] 




D 





0(2) 



D = p 2 (l + [rn f + a s )e 



aL;. 



+ z- 



'([m f + a s 



-aL. 



) + \n s \ 2 (z 2 +p 2 e 



2aL. 



If L s is an even number or if m is such that < b(p) (and therefore < e a ^ 

-aL 



(34) 

(35) 

(36) 
) for 

all momenta then e~ aLs is positive and my must be negative in order to have a phase 
with 7r s 7^ 0. The phase boundary mj^ft) can be obtained from the above equations by 
setting 7r s — > and eliminating x = mj + a s . From these equations it can be deduced 
that \irif c \ decreases exponentially with L s . Furthermore, for a given (3 lc the width of the 
7r s 7^ region also decreases exponentially with L s . In figure 4 the phase boundary of the 
7r s phase is given for a 6 3 lattice with antiperiodic boundary conditions along the time 
direction, m = 1 and for various L s values. From bottom to top the L s values are 2, 3, 4, 
and 5. This has been calculated using the full form of the SPE and not just the truncated 
form eq. 



34]and|35|. However, the two are nearly the same for 2 < L s . Also different values 
of mo produce similar results. As can be seen from figure 4 when L s is increased the phase 
boundary moves to smaller \rrif\ with decreasing width. In figure 5 for fixed j3\ = 0.05 the 
origin of the phase boundary —rrif c and width 5m f c are plotted versus L s for mo = 1.0 and 
a 6 3 lattice with antiperiodic boundary conditions along the time direction. 

The effective fermion mass m q is identified from the zero of eq. for momenta p ~ 
and is: 



m n 



z(p = 0)[m f + a s + (1 - m ) 



(37) 



The continuum limit is reached at m q = and the lattice spacing is set to a ~ m q . This 
point corresponds to the largest /3i c value of the phase boundary curves in figure 4 [m q 
is positive on the upper part of the phase boundary and negative on the lower one). The 
width of the Aoki phase close to m q = can be obtained from eq. [J3| To lowest order in m q 
is: 5mj = m q ^f\ mq= o. The derivative is not zero at m q = and therefore 5mj ~ m q ~ a. 
For an analysis concerning the width of the Aoki phase in QCD with Wilson fermions see 
27| . The above features are not particular to three dimensions and similar results have 



been obtained for four- dimensional lattices. 

Finally, it is interesting to observe that if m is such that b(p) < for some momenta 
and L s is odd then the 7r s 7^ phase can occur even for positive mj. If the 7r s 7^ phase 
needs to be avoided for any mo then one should set < rrif and L s to an even number (on 
the other hand, if mo < 1 then < b(p) for any momenta and any value of L s can be used). 



4 Hybrid Monte Carlo simulations 

In this section standard Hybrid Monte Carlo (HMC) simulations are performed for N = 2. 
These simulations support the large iV results of the previous section. For all simulations 
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the trajectory length is set to r = 1 and the step size to 5t = 0.1. The acceptance rate is 
~ 90% and the conjugate gradient inverter iterations are ~ 50 — 100. Typically 100 — 200 
thermalization sweeps were followed by 300 — 400 sweeps with measurements. The lattice 
size for all simulations is 6 3 with antiperiodic boundary conditions along the time direction. 
All simulations were done on workstations. 

The Z2 x Z2 model with action as in eq. |l] and interaction term as in eq. 13 was 
simulated first. The initial configuration for all simulations was a uniform configuration 
with a — 1. In figure 6, < a > is plotted as a function of /3i for mj = 0.02, L s = 12, 
m = 0.4 (diamonds), m = 1.0 (squares) and m = 1.6 (crosses). The solid lines are the 
large N predictions for the same parameters and the dotted lines are the large N predictions 
for rrif = and L s = 00. The agreement with the large N predictions is good away from 
the critical region where the approximation of a uniform saddle is good. 

The decay and decay rate of < o > as a function of L s for rrif = 0.02 is shown in figure 
7. In figure 7a m = 1.0 and (3\ = 0.05 (diamonds), fti = 0.1 (squares) and /3i = 0.25 
(crosses). In figure 7c (3i = 0.05 and m = 0.6 (diamonds), mo = 0.8 (squares), m = 1.0 
(crosses), and mo = 1.1 (circles). The solid lines in figures 7a and 7c are fits to co + cie _C2is 
and the dotted lines are the large N predictions. The decay rates C2 from the fits in figures 
7a and 7c are shown in figures 7b and 7d respectively. The solid lines in these figures are 
the large N predicted minimum and maximum decay rates. As can be seen, the agreement 
with the large iV predictions is good and the decay rates are fairly independent of (3\ and 
mo and close to the minimum predicted value. This is in contrast to gauge theories where 
the dependence of the decay rates on the gauge coupling is significant 0, [| . 

Finally the SU(2) x 577(2) model with action as in eq. |I] and interaction term as in 
eq. |6| was simulated in an effort to investigate the presence of the parity-flavor broken 
phase. The lattice Lagrangian has exact SU (2) flavor symmetry and as a result when this 
symmetry is spontaneously broken there will be two exactly massless Goldstone pions. In 
that case the Conjugate Gradient inverter would not converge. Furthermore, in the small 
lattices considered here spontaneous breaking can not really take place and the tz field 
would always average to zero for sufficiently large statistics. In f22"l , |28| , f29f these problems 
were treated by adding a small external field h that breaks the SU(2) flavor symmetry and 
therefore not only gives mass to the two pions but also provides a "preferred" orientation 
for the 7r field. Such an external field is used here by adding a term that is exactly the 



same as the 7r term in eq. 16 but with n replaced by h. The initial configuration for 



all simulations was a uniform configuration with (a, tti, tt 2 , vr 3 ) = (1,0,0,1). In figure 8 
the average value of 7r 3 is plotted versus rrif for a 6 3 lattice with antiperiodic boundary 
conditions along the time direction, h = (0,0,0.1), j3\ = 0.05, L s = 2 and m = 1.0. 
The "outer" solid line is the large N prediction for h = (0,0,0.1) and the "inner" one for 
h = (0,0,0). The diamonds are the results of the numerical simulations. As can be seen 
they are in fairly good agreement with the large N results supporting the presence of this 
phase. A finite volume analysis together with an analysis involving decreasing values of h 
as in |29j is still needed in order to unequivocally establish the presence of the parity-flavor 
broken phase for N = 2, but this is not in the scope of this paper. 
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5 Conclusions 



The interaction of domain wall fermions with scalar fields was formulated. It was found 
that contrary to naive expectations this interaction takes place only at the boundary of the 
extra direction. This is in contrast to the interaction of domain wall fermions with gauge 
fields which is the same along the extra direction. This seems to indicate a picture with 
some richness where different spin fields couple to domain wall fermions in different ways. 

Large N techniques were used to study two 4-Fermi models, one with Z 2 x Z 2 and 
one with SU(2) x SU (2) chiral symmetry. It was found that at the limit of infinite extra 
direction the chiral symmetry breaks spontaneously in the standard fashion. However, if 
the size of the extra direction is finite the SU(2) x SU(2) chiral symmetry is explicitly 
broken by the regulator down to flavor SU(2). It was found that this remaining flavor 
symmetry can break spontaneously if the bare quark mass is negative resulting to a parity- 
flavor broken phase of the Aoki type. Hybrid Monte Carlo simulations were performed for 
those models with iV = 2 on 6 3 lattices with antiperiodic boundary conditions along the 
extra direction. The results were found to support the large N predictions. 
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Figure 1. The large N critical value lc versus m for a 6 3 lattice with antiperiodic 
boundary conditions along the time direction. 




0.0 0.5 1.0 1.5 2.0 



m 

Figure 2. The real part of the maximum and minimum large N decay rates as a function 
of m for a 6 3 lattice with antiperiodic boundary conditions along the time direction. 
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Figure 3. The large N a s as a function of L s for a 6 3 lattice with antiperiodic boundary 
conditions along the time direction, m,f — and various m at (5\ = 0.3 above the transition. 
The diamonds, squares, crosses, circles, and stars correspond to m = 0.2,0.4,0.6,0.8, 1.0. 




Figure 4. The large N phase boundary of the parity-flavor broken phase of the SU (2) x 
SU (2) model on a 6 3 lattice with antiperiodic boundary conditions along the time direction 
with itlq = 1 and for various L s values. From bottom to top the L s values are 2, 3, 4, and 
5. The parity-flavor symmetry is broken inside the oval-looking regions. 
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Figure 5. For fixed 0i = 0.05 the origin -mj c (diamonds) and width 5m f c (circles) of 
the large N parity-flavor broken phase of the SU (2) x SU (2) model is plotted versus L s for 
m = 1 and for a 6 3 lattice with antiperiodic boundary conditions along the time direction. 
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Figure 6. < a > as a function of f3\ for the Z2 x Z2 model on a 6 3 lattice with antiperiodic 
boundary conditions along the time direction, mj = 0.02, L s = 12, m = 0.4 (diamonds), 
m = 1.0 (squares) and m = 1.6 (crosses). The solid lines are the large N predictions 
for the same parameters and the dotted lines are the large N predictions for mj = and 
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Figure 7. The decay and decay rate of < a > as a function of L s for the Z 2 x Z 2 model on 
a 6 3 lattice with antiperiodic boundary conditions along the time direction and nif — 0.02. 
In figure 7a m Q = 1.0 and j3i = 0.05 (diamonds), (3\ — 0.1 (squares) and j3i = 0.25 (crosses). 
In figure 7c /3i = 0.05 and m Q = 0.6 (diamonds), m = 0.8 (squares), m = 1.0 (crosses), 
and m = 1.1 (circles). The solid lines in figures 7a and 7c are fits to c + C\e~ C2hs and the 
dotted lines are the large N predictions. The decay rates c 2 from the fits in figures 7a and 
7c are shown in figures 7b and 7d respectively. The solid lines in these figures are the large 
N predicted minimum and maximum decay rates. 
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Figure 8. The average value of the third component of the pion field versus mf for the 
577(2) x SU(2) model on a 6 3 lattice with antiperiodic boundary conditions along the time 
direction, h = (0,0,0.1), fa = 0.05, L s = 2 and m = 1.0. The "outer" solid line is the 
large N prediction for h = (0,0,0.1) and the "inner" one for h = (0,0,0). The diamonds 
are the results of the numerical simulations. 
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